Anomalous decay of a prepared state due to non-Ohmic coupling to the continuum 
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We study the decay of a prepared state Eo into a continuum {Ek} in the case of non-Olimic 
models. This means that the coupling is |Vfc_o| oc \Ek—Eo\''~^ with s 7^ 1. We find that irrespective 
of model details there is a universal generalized Wigner time to that characterizes the evolution 
of the survival probability Po{t)- The generic decay behavior which is implied by rate equation 
phenomenology is a slowing down stretched exponential, reflecting the gradual resolution of the 
bandprofile. But depending on non-universal features of the model a power-law decay might take 
over: it is only for an Ohmic coupling to the continuum that we get a robust exponential decay 
that is insensitive to the nature of the intra-continuum couplings. The analysis highlights the co- 
existence of perturbative and non-perturbative features in the dynamics. It turns out that there are 
special circumstances in which to is reflected in the spreading process and not only in the survival 
probability, contrary to the naive linear response theory expectation. 
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I. INTRODUCTION 

The time relaxation of a quantum-mechanical prepared 
state into a continuum due to some residual interaction 
is of great interest in many fields of physics. Applications 
can be found in areas as diverse as nuclear [l| , atomic and 
molecular physics_0 to quantum information 0], solid- 
state physics 0, 01 and quantum chaos Q. The most 
fundamental measure characterizing the time relaxation 
process is the so-called survival probability Po(t), defined 
as the probability not to decay before time t. 

The study of Po{t) goes back to the work of Wcisskopf 
and Wigner 0, [1] regarding the decay of a bound state 
into a continuum. They have found that Po{t) follows 
an exponential decay Po{t) = exp(— i/io), with rate l/tg 
given by the Fermi Golden Rule (FGR). 

Following Wigner, many studies have adopted Ran- 
dom Matrix Theory (RMT) modeling for the in- 
vestigation of Poit), highlighting the iinportance of the 
statistical properties of the spectrum both in the 
semiclassical ^^'^ i^. the many-body context [Tsl 
where non-exponential decay such as Po{t) ^ cxp{—^/i) 
may arise. 

Despite the interest in specific problems where devia- 
tions from the Wigner theory arise, a theoretical inves- 
tigation of the time relaxation for prototype RMT mod- 
els is still missing, and also the general (not model spe- 
cific) perspective is lacking. This unbalanced situation 
should be contrasted with the arena of spectral statis- 
tics [l^, where one can find on the one hand elaborated 
mathematical studies of generics RMT models, and on 
the other hand system specific investigations that explore 
non-universal features that cannot be captured by RMT. 

In this paper we would like to explore the limits of 
universality in decay problems using an RMT perspec- 
tive. Specifically we explore the decay of an initially 
prepared state into the continuum assuming non-Ohmic 
rather than Ohmic circumstances. The latter notions 
are precisely defined in the next section. We show that 
the survival probability Po{t) = g{t/to) is characterized 



by a generalized Wigner decay time Iq that depends in 
a non-linear way on the strength of the coupling. We 
also establish that the scaling function g has distinct uni- 
versal and non-universal features. It is only for Ohmic 
coupling to the continuum that we get a robust exponen- 
tial decay that is insensitive to the nature of the intra- 
continuum couplings. In addition to Po{t) we investigate 
other characteristics of the evolving wavepacket: the vari- 
ance A£',p,d(i) and the 50% probability width AE,„,,{t) 
of the energy distribution, that describe universal and 
non- universal features of its decaying component. 

The structure of the paper is as follows: In Section II 
wc define the Friedrichs model (FM) [20| and the gener- 
alized Wigner model (WM), and discuss the numerical 
approach used in the subsequent sections. The quanti- 
ties under investigation, and their physical meaning is 
discussed in Section III. In Section IV we present an 
overview of the relevant time scales that dictate the dy- 
namics of our models. In Section V we present analyti- 
cal and numerical results for the Local Density of States 
(LDOS) of the FM and of the WM. The study of the 
wavepacket dynamics is presented in Sections VTVIII: 
first we analyze the decay of the survival probability 
and later the evolution of the energy distributioir. Our 
conclusions are given at the last section IX, where we 
also discuss the crossover from the universal to the non- 
universal behavior. 



II. MODELING 

We analyze two models whose dynamics is generated 
by an Hamiltonian 



V 



(1) 



with Hq = diagji^n} and n G Z. The system is prepared 
initially in the eigenstate corresponding to Eq, and the 
coupling to the other levels is characterized by the spec- 



2 



tral function 

CH = ^|K,o|'2^,5(c^-(i;„-£;o)) (2) 

In the so-called Ohmic case C{uj) is finite for w ^ and 
accordingly FGR suggests a well definite finite rate of 
decay. But our interest is in so-called non-Ohmic cir- 
cumstances for which C{uj) either vanishes or diverges as 
w — >■ 0. To be specific we assume the standard form 



C(w) = 27re2|w| 



(3) 



Note that C{oj) has 1/time units, and hence e has 
f/time^~* units. In order to define the model one should 
specify the exponent s, the strength of the coupling e, 
the density of states g and the bandwidth = hQ~^. 
But this is not enough. One also has to specify the cou- 
plings between the other levels En-tQ. Here we distinguish 
between two cases that arc illustrated in Fig.(IT]) and dis- 
cussed below. 

The Friedrichs model (FM) features a distin- 

guished energy level Eq that is coupled to the rest of the 
levels En^Q by a rank two matrix. This means that the 
other levels are not coupled. For example the following 
is a 6 X 6 FM matrix with 6 = 2. 



/ Eo Vol Vo2 
Vio -El 



n 



V20 





V 



E2 
E3 
E4 



\ 









E5 / 



In the FM case the dimensionless bandwidth b defines the 
effective size of the matrix as = 6 -I- 1, because out-of- 
band levels cannot be reached. 

The Wigner model (WM) 0] features a perturbation 
matrix that does not discriminate between the levels, and 
is given by a banded random matrix. For example: 



I Eo Vol V02 \ 

Vio El V12 Vi3 

V20 V21 E2 V23 V24 

1^31 V32 E3 Vu V35 

V32 V43 E4 V45 

V 1/43 144 E5 1 



(5) 



Here we assume that all the levels arc alike, and therefore 
Eq.([3]) implies that 



IK 



oc 



\Er,~Er, 



(6) 



where the overline indicates an RMT averaging over real- 
izations, which from now on will be implicit in the analy- 
sis of the WM case. Unlike the FM case, here the dimen- 
sion of the matrix N has significance, because any level 
can be reached by high order processes. As far as the 
theoretical analysis is concerned we assume the matrix 
to be of infinite size. 



The assumed form Eq.([3|) for the spectral function 
C{bS) is motivated by the study of various model systems 
(e.g. those of Refs. |l3 - [l6l [l8l| ) where the assumption of 
"flat" band-profile looks like an oversimplification. Thus 
it constitutes the natural generalization for the standard 
FM and WM. By integrating Eq.© over w we see that 
the perturbation V that appears in Eq.© is bounded 
provided s > 0. The s = 1 case is what we refer to as the 
Ohmic case, for which it is well known that both mod- 
els leads to the same exponential decay for the survival 
probability Q . For s > 2 the effect of the continuum can 
be handled using 1°' order perturbation theory. We focus 
in the < s < 2 regime and consider the s ^ 1 case for 
which a non-linear version of the Wigner decay problem 
is encountered. 

We measure the energy taking Eo reference 
and accordingly we use the notations lo = E — Eq, and 
ujn = E„ ~ Eo- In the continuum limit any summation 
over uj„ is replaced by an integral over a variable uj' . See 
Fig.® for illustration. 

In the numerical simulations we integrate the 
Schrodinger equation for the amplitudes ipn{t) = ('T'IV'(O) 
starting with the initial condition V'n(O) = <^ti.o- We use 
units such that h=l, the density of states is g~l, and 
Eq—O, and we assume a sharp bandwidth 



b = gujc 



[bandwidth] 



(7) 



The matrix elements of V are taken from a Gaussian dis- 
(4) tribution with zero mean and variance in the units 
as defined above. The integration is done using the self- 
expanding algorithm of 21| to eliminate finite-size effects, 
adding 106 sites to each edge of the energy lattice when- 
ever the probability of finding the 'particle' at the edge 
sites exceeds 10~^^. 



III. STRATEGY OF ANALYSIS 

We denote by |n) the unperturbed eigenstates of Ho 
and by \E^) the perturbed eigenstates of H. The local 
density of states LDOS with respect to the initial state 
|0) is defined as 



(8) 



Starting with the initial state the time-dependent state 
lipit)) can be represented by the amplitudes c„(t), such 
that 



e-'-^"* In) 



The probability distribution at time t is 

P„(i) = \{n\Uim\^ = |c„(i)P 



(9) 



(10) 



where U{t) = cxp[—iHt] is the evolution operator. The 
survival probability is Po{t). We can associate with the 
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energy distribution a probability density 

PtH = ^P„(i) ,5(w-c^„) (11) 

n 

The energy distribution is characterized by its dispersion 

-I 1/2 



(12) 



and by the median £^50% = Eq, and also by the i?25% 
and £'75% percentiles. The width of the core component 
is defined as 

^Ecoxc{t) = £'75% - i?25% (13) 

Our interest below is focused in Po(i)i in A£'^o„(t) and in 
A£^Bprd(0- We use the same measures in order to describe 
the LDOS distribution p{ijj). See Fig. ([3]) for a cartoon 
that illustrates the significance of the different measures 
in the analysis. 

The first stage of the analysis is to find the LDOS, 
also known as the strength function. In the FM case it 
can be done analytically using a standard Green function 
method. In the WM case it is possible to generalize the 
approach by Wigner and followers [1, Q . 

The survival amplitude is related to the LDOS. 
Namely, Co(i) can be written as the Fourier transform 
(FT) of p{u)). The derivation of this well known fact is 
simple. Taking the energy reference as E'q = wc have: 



co(<) = (0|;7(t)|0) 



-mt 



E,){E„\Q) 



,2 ^-^E^t 



FT[27rp(a;) 



(14) 

(15) 
(16) 
(17) 



For obvious reasons there is resemblance between the 
saturation profile Poo('^) and the LDOS p(w). In a 
stochastic (diagonal) approximation the former is the 
auto-convolution of the latter. But in order to find the 
time dependence of AE^aio{t) and AE^-p^ait) it is not 
enough to know the LDOS. In order to obtain analytical 
results for the spreading we shall use a linear-response 
strategy that can be further refined in the FM case. 



IV. TIME SCALES 

There are two frequency cut-offs, an infrared cutoff 
LOg ~ g^^ and an ultraviolet cutoff which are the mean 
level spacing and the range of the coupling, respectively. 
The associated time scales are the Heisenberg time tn 
and the semi-classical time tc which are given by 



tc 



2ttq 



(18) 
(19) 



We shall see that the continuum limit p"-"^ — )■ is well 
defined if s > 0. If we further assume s < 2 then well 
defined results are obtained also in the universal limit 
— > 00. Thus in the range < s < 2 we should have in 
the continuum limit a cutoff free universal theory, that 
constitutes a generalization of the Wigner decay problem. 
Note that for finite values of the cutoffs the actual range 
is more restricted, as explained in Sect IV Al and depicted 
in Fig.®. 

We shall see that the decay is characterized by what 
we call generalized Wigner time: 



tn = 



2^e2 



r(3-s)sin(s7r/2) 



-1/(2- 



(20) 



where T is the Gamma function. The numerical prefac- 
tor is explained and calculated in Sect. (| VI Ap . As far as 
order-of-magnitude estimates are concerned a more prac- 
tical expression is 



2^,2 X -1/(2-.) 



(2-s)s 



(21) 



One observes clearly that s— 7>0 and s— 7>2 are limiting cases 
that require special attention due to the increased sensi- 
tivity to the infrared and to the ultraviolate cutoffs re- 
spectively. 



V. THE LOCAL DENSITY OF STATES 

Following the presentation of [13, HH we expect the 
LDOS to have in the most general case three regions: 
(i) The core that consists of the levels that are mixed 
non-perturbatively; (ii) The first order tails; (iii) The 
higher order tails. We now explain this terminology and 
the associated phenomenology. 

The LDOS is determined by the overlaps {E,y\0) of 
the initial state with the perturbed eigenstates. So nat- 
urally the first question that arises is whether we can 
use perturbation theory for the calculation. First order 
perturbation theory (FOPT) assumes that we can as- 
sociate perturbed and unperturbed levels (no mixing), 
hence {E^=Q\n=0) ~ 0(1). If this is indeed the case we 
say that the core consists of one level only, while the tails 
consist of all the {E^^q\ji—0) overlaps. Typically FOPT 
gives a leading order approximation for the tails, while 
higher orders are essential for those levels that are not 
coupled directly [T^] . In our model the higher order tails 
reside at uj > ujc out of the range of physical interest, 
because lOc is assumed to be very large. 

The first order tails are fully determined by the spec- 
tral function C{uj). The self consistent condition for the 
validity of FOPT is po ^ 1 where po is defined as the 
probability which is carried by the tails: 



En — Eq 



C{uj) dw 
uj'^ 2^ 



(22) 
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By substituting the spectral function of Eq.(l3l) we get calculation (App.B). The final result is 



Po 



2 /2r(s-2)w/-2 foj. 



(23) 



In the continuum limit po is infrared divergent if s < 2, 
and consequently FOPT does not apply even if e is made 
very small. This is the case of our interest, and it should 
be contrasted with the s > 2 case for which FOPT applies 
in the weak coupling limit. In the latter case the LDOS 
can be written schematically as 



p{uj) {l-po)S{u) + PQ—f( — 



(24) 



where the second term stands for C (u) / {2ttuj'^) . The 
natural question that arises is what happens to this line 
shape in the s < 2 case, and in particular whether there 
are remnants of FOPT. It turns out that generically the 
answer is positive [23| . In spite of the mixing of nearby 
levels, the tails are still given by the FOPT expression. 
Accordingly we can identify in the LDOS a core region 
1^1 < 7o that contains the large overlaps, and FOPT tails 
that dominate the outer |a;| > 70 regions. The char- 
acteristic frequency can be determined either analyti- 
cally or self-consistently and consequently the generalized 
Wigner time is identified as to = I/70. This phenomeno- 
logical picture will be formulated in a more precise way 
in the subsequent subsections. 

We are going to analyze the LDOS for s < 2, where 
we have to go beyond FOPT, but still can use the core- 
tail phenomenology. The non-perturbative core within 
\u)\ < 7o has non-universal structure that depends on the 
details of the model. There are various strategies to de- 
duce the structure of the core. In particular: the Green 
function method that we are going to use in the FM case; 
the RMT calculation of the moments as in the pioneering 
work of Wigner [1, Q ; and the semiclassical reasoning if 
the model has a classical limit 

A few words about the Lamb shift are in order. FOPT 
allows to calculate the shift of Eq due to the repulsion by 
the other levels: 



A(0) = E 



\Vn,0\ 



Eq — En 



C{uj') duo' 
u' '2tt 



(25) 



If s > the lamb shift is not infrared divergent, provided 
we keep away from the energy floor, and we get A(0) = 
due to the symmetry of (^(a;). This of course does not 
mean that Lamb-shift physics is irrelevant. In the anal- 
ysis of the FM model we shall define a spectral function 
A(aj) that plays a major role in the analysis. 



A. The LDOS - FM 

In the FM case it is possible to derive an exact ex- 
pression for the LDOS either via the Green function for- 
malism (App.A) or from a straightforward elementary 



where 



r(c.) 

A(a.) = 



r(^)/2 



IT (C^ - A(t^))2 + (r(c^)/2)2 



E|K,o|'27r,5(S-S„) = (7(c 



IK 



n,0| 



E — En 



+ 00 



C{ijj') div' 
oj — ui' 27r 



(26) 

(27) 
(28) 



In the Ohmic case (s=l) we have A(cj) = 0, and 
r(a;) = const, and consequently the LDOS is the familiar 
Lorentzian. In order to calculate A.{lu) for s 7^ 1 we ex- 
ploit the fact that C{uj) is an even function, while A(aj) 
comes out odd. Gonsequently we can write Ea.p8| as 



A{uj) = - 

TT 



{uj'y-^doj' 



(29) 



-e^jwl'* -^sgn(w) 



ln|tJc/tj| (s-l)a: 



dx 



where 



is the level spacing, and we used the 



substitution uj' = cje^. For 0<s<2 it is possible to take 

00. Using the integral 



the limits cjg— >0 and oj, 



+00 



sinli((s — l)x) 



dx 



-7rcot(s7r/2) (30) 



, sinh(a;) 

we get the cutoff-free result 

A{u) = e^TT cot (s7r/2) |w|'"^sgn(w) (31) 

In contrast to that, the marginal cases 5=0 and s—2 are 
lower and upper cut-off dependent respectively: 



A(a.) 
A(a.) 



-e^ LU In 



1 In 



1 - {uJc/ujY 



(Uj/Ujg) 



1 



s=2 
s=0 



(32) 
(33) 



The crossover from the "0<s<l" to the "s=0" and the 
"s=2" expressions is not sharp. Looking carefully at the 
integral we see that the condition for a cutoff independent 
result is (see FigH]): 



(34) 



Having obtained A(a;) we can substitute it into Eq. ipS)) 
and get the LDOS. Both F and A are cx u>^~^, with 
s dependent prefactors. One can define a characteris- 
tic crossover frequency 70 above which the term in 
the denominator if Ea. (j26|) dominates. In the universal 
regime this leads to 



70 



sin(s7r/2)| 



l/{2-8) 



(35) 
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while in the hmiting cases it attains the values 

70 ~ Wee" as s^2 (36) 
70 ~ eln{ge) as s— >0 (37) 

Depending whether uj ^ or w 3> 7o one obtains dif- 
ferent approximations for the LDOS that can be packed 
together in the following writing style 



1 



r/2 



r/2 



(r/2)^ 



Consequently we get for 0<s<2 



'sin(s7r/2) 



(38) 



(39) 



This validity regions of the expression in the curly brack- 
ets are illustrated in Fig.(|3]). After FT they are reflected 
in time domain as implied by the analysis in Fig.([5|), 
which we further discuss later. The numerically deter- 
mined LDOS in the FM case is displayed in Fig.®, and 
contrasted with the WM case. 

Eq.([39|) holds provided 70 is in the range that has been 
specified in Eq. (j34[) . This constitutes a condition on e 
that can be regarded as a generalized version of the FGR 
condition. For marginal values of s it becomes more diffi- 
cult to satisfy the generalized FGR condition, and even- 
tually the range of the universal behavior as defined in 
Ea. (p4)) shrinks to zero (see Fig.Q). In particular for 
s— >2, once the generalized FGR condition breaks down, 
we get 



Pi-) - 73 



1 



1 



4e2 \uj\ln^ \oj,/u;\' \oj\ 



(40) 



In this expression the l/co tails prevail for a; > 70, anal- 
ogous to Eg. ([5^ . but here 70 is ujc dependent. 



B. The LDOS - WM 

As discussed in the beginning of this section, for s > 2 
we can calculate the LDOS using FOPT. In contrast, 
for s < 2, FOPT breaks down for any finite coupling, 
because in the continuum limit po is infrared divergent. 
Still we can use the FOPT result down to a frequency 
7o which we estimate below in a self consistent man- 
ner. For lower frequencies the LDOS is dominated by a 
non-perturbative core which is model specific and in that 
sense, non-universal. For the WM we write schematically 



Semi-circle, 



|3-s 



(41) 



This functional form is illustrated in Fig. ([3]), and tested 
numerically in Fig.®. The reasoning behind this expres- 
sion is further detailed below. 



The crossover frequency 70 is almost the same as 70 
which we define in the time domain analysis of the next 
section [the distinction between the two is clarified after 
Eg. ([55)) ]. We define the crossover frequency 70 by the 
equation 



which gives 



|i^|>7o 



7o 



C(a;) du! 
2^ 



50% 



2 \ l/(2-s) 



2-s 



(42) 



(43) 



The FOPT tails are dominated by "Ho, while V is re- 
garded as a perturbation. In contrast to that the core is 
fully dominated by V . Following Ref.[^ we can argued 
that the core of the LDOS is a semicircle of width AE^^, 
which in our case is given by the expression 



1/2 



(44) 



In this expression we use the effective bandwidth 70 
rather than the real bandwidth ujc, as implied by our 
core-tail hypothesis. Performing the integration we get 



7o 



(45) 



which implies that the calculation is indeed self- 
consistent, i.e. there is a well defined frequency 70 at 
which the core-tail crossover takes place. 



VI. THE DECAY OF THE SURVIVAL 
PROBABILITY 

We turn now to the study of wavepacket dynamics, 
with the objective to contrast the results of the FM with 
those of the generalized WM. We assume that the initial 
preparation is delta-like in energy space. The tempo- 
ral behavior of the wavepacket, is characterized by three 
quantities: (a) the survival probability is directly related 
to the LDOS analysis; (b) the core width refiects the 
non- universal component of the wavepacket; (c) the en- 
ergy spreading refiects the tails of the wavepacket. 

In the present section we study the survival ampli- 
tude, which is obtained via FT of the LDOS, namely 
co(t) = FT[27rp(a;)]. For s>2 the FT of the FOPT esti- 
mate Eg. ((M)) has a simple crossover at tc to saturation 



co(t) 



(1-Po) 
1 -Po 



^P0f{t/tc) 
1- ' 



(46) 



For s<2 the core is not a discrete delta function but 
rather has a width 70 which implies a progressive de- 
cay on time scale l/7o- Our focus below is on the regime 
0<s<2, and we pay some extra attention to the limiting 
cases s— and s— )■ 2. 
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A. The survival probability - WM 



expression for the crossover time 



In the WM case the LDOS has power-law tails 
l/jwl^"''" with a = 2— s, and a smooth semicircle core. 
Therefore its FT is exponential-like for all times. The 
similarity to the a-stable Levy distribution suggests a 
stretched exponential decay which is confirmed by our 
simulations. There is an optional argument that sup- 
ports a stretched exponential behaviour. The stretched 
exponential can be regarded as the solution of the fol- 
lowing phenomenological rate equation for the survival 
probability 



V{t) = -r(const/i) Po{t) 



(47) 



This phenomenological equation involves a time depen- 
dent effective decay rate into the quasi continuum. It 
is estimated as r(a;) with w ~ 1/t which reflects the en- 
ergy uncertainty at time t. This phenomenological ap- 
proach does not provide a definite result for the numerical 
const. A reasonable procedure is to determine this con- 
stant from the short time behavior, which is determined 
by the tails. We recall the following FT relation: 



|wr cos[ujt) — 

2TT 



ir(/3) cos(/3| 



1 



where /3 > 0. Switching the roles of t and w, the inverse 
FT implies that power law tails l/lwl^"*"" are related to 
discontinuity or singularity C\t\" at time domain, where 
C = -[2r(l-ha)sin(a7r/2)]-^ This leads to the resuh 



co{t) 



exp 



2-s 



(49) 



where to is defined as in Eq. (pO)) . Going back to the phe- 
nomenological rate equation it implies that the numerical 
constant is [(2-s)/(r(3-s) sin(s7r/2))]i/("-i). 



B. The survival probability - FM 

In the FM case the core is singular and therefore a 
power-law takes over for t > too, where the crossover time 
too is determined below. In Fig. ([5]) we present numerical 
FTs that support this statement. Note that in the s=2 
limit the power law becomes a logarithmic decay. 

In order to obtain the explicit expression for the power 
law decay we use the FT relation of Eq. p8)) . Due to 
the discontinuity or singularity at w = we get for large 
times 



coit) 



sin((s-l)7r) sin(s7r/2)r(2- 

2sin((s-l)7r) Ao " 
{2-s)n \T 



(50) 



where in the last equation we expressed co(t) by tg. Com- 
paring the exponential-like and the power law we find an 



2 In 



2sin(|s-l|7r) 

{2-s)n 



l/(2-s) 



in 



(51) 



This expression diverges in the limit s = 1 , implying that 
only the exponential survives. The numerical simulations 
presented in Figs.© and ^ support our findings. 

The marginal case s=2 requires special treatment. In 
this case the strength of the interaction (e) is a dimen- 
sionless parameter. The oc l/w tail of the LDOS in 
Ea. (|40p is ujc independent, but with the uJc dependent 
lower cutoff of Eq. ([55)) . For the purpose of FT we ap- 
proximate the LDOS as 



p{uj) 



(52) 



and from the integral representation of the modified 
Bessel function of the second kind, using the approxi- 
mation Ko{t) « const — ln(T) it follows that 



co{t) 



log 



for < i < to (53) 



In order to clarify how the s=2 case is related to the gen- 
eral s expression, it is useful to note that for very small /3 
one has t^ ~ /3\n{t). It is also important to realize that 
for t < to both the s=l exponential and the s=2 log- 
arithm are consistent with the short time 1 — (t/to)^~^ 
expression that holds for general s. In contrast to that 
the long time decay is sensitive to the small frequency 
features of the LDOS. 



VII. THE EVOLUTION OF THE CORE WIDTH 

In this section we study the time evolution of the core 
as refiected in the width of the energy distribution. The 
width AE^^^^{t) is intimately related to Po{t). It starts to 
rise at t ~ when Po{t) = 50%, and its saturation value 
is determined by the core width of the LDOS ^ 70. Thus 
AE^„^ should exhibit one parameter scaling with respect 
to to. This has been verified numerically, and is shown 
in Fig.®. A more elaborated analysis follows below. 

For i <C to the initial state carries most of the probabil- 
ity, and therefore we can use FOPT to obtain in analogy 
with Eq. ((M)) : 



Pt(w) 



Sic.) 



1 (7(0;) 
2^ w2 



[2sin(wt/2)]2 (54) 



with normalization measure dio. As long as this expres- 
sion holds we say that the core component contains one 
level only. It is clear that after a longer time pt{uj) be- 
comes similar to the LDOS p(o;). 

Accordingly it is natural to ask what is the time evo- 
lution of the core-tail structure. In particular one would 
like to estimate the time dependence of the border 7(t) 



7 



between the core and the taiL Using the same reasoning 
as in Eq. p^ we get for 7(t) the equation 



s 2— s 



50% (55) 



From this equation it foUows that there is a sharp 
crossover from -f{t) = g^^ to '-f{t) ^ 7o during the time 
interval [to, to] where tg is given by Ea. (pT|) and to = I/70 
is given by Eq. p3)) . In Fig.® we present the resuhs of 
the numerical analysis. Our data, indicate that the ex- 
pected one-parameter scaling is obeyed. We have verified 
that the deviations from the expected e dependence (for 
large e) are an artifact due to having finite (rather then 
infinite) bandwidth. 



VIII. THE TIME DEPENDENT SPREADING 

In this section we analyze a very different character- 
istic of the evolving wavepacket, which is the spreading 
A£'3p,d(i). Unlike the width which is determined by the 
core, the spreading is determined by the tails of the distri- 
bution. Irrespective of whether we deal with FM or with 
WM the tails are determined by FOPT and therefore the 
naive expectation is to have a tc dependent rather than 
to dependent evolution. We shall see that this is roughly 
but not quite correct. 

For the analysis one can use the traditional strategy 
of Refs.[l3, Hll, [H], leading to the following spreading 
formula (see App.C): 

[^E^,Atf = C(0, 0) + Cit, t) - C(0, t) - Cit, 0) (56) 



where 



C{t',t") = {V{t')V{t")) 



(57) 



This expression is formally exact. But in order to get 
a practical expression the LRT approximation assumes 
that it is possible to calculate the correlation func- 
tion with the unperturbed Hamiltonian. Accordingly 
C{t',t") C{t' -t") where C{t) is the inverse FT of 
C{oj), leading to 



!(C(0)-C(i) 



This expression implies that 
AS,,,, (00) 



2 X 

s 



1/2 



1/2 



(58) 



(59) 



This saturation values, which diverges in the Wc — ^ c» 
limit, is attained after a short transient time tc ^ 0, 
while the generalized Wigner time to is not reflected. The 
above expressions agrees with the numerical calculations, 
as shown in Fig. pHj) . 

The LRT approximation C{t\t") w C{t' - t") makes 
sense if the correlation function is stationary. This is not 



the case with FM. By inspection a more appropriate ap- 
proximation in the latter case should take into account 
the decay of the initial state and hence it should incor- 
porate the survival amplitude co(t) as follows: 

C{t',t") « co{t')co{t")C{t' -t") (60) 

Within the framework of this approximation we get 

1/2 



AS=p,d(t) « (l+Po(i))C(0) - 2co{t)C{t) 



(61) 



where Po{t) = |co(t)p. This leads to a saturation value 
smaller by factor compared with the WM case, re- 
flecting the non-stationary decay of the fluctuations as a 
function of time. More interestingly Ea. (j6ip contains a 
cutoff independent term that reflects the universal time 
scale to- 

In fact in the FM case it is possible to obtain an exact 
result. First we notice that by deflnition co(t) — (e~'^*), 
where the expectation value here and below is with re- 
spect to the initial |0) state. By differentiating we get 
either co(t) = -i('He-''^*), or co(t) = -i{c-'''^^ H) . Set- 
ting Ti, = Ho + V and using the convention Eq—Q for the 
|0) eigenstate of "Hq we find 

co{t) = -i{Vc-'^') = -i{c-'^'V) (62) 

Similarly we can handle the second derivative leading to 

co(t) = -{Hc-'^'H) = -{Vc-''^*V) (63) 

We now realize that the correlations in Eg. ([55)1 can be 
expressed using C{t) and co(t). Trivially we have 



{vmm = {V') = c(o) 



(64) 



In the FM model the only nonzero elements of V are Vb,n 
and Vnfl with n ^ 0. Consequently we can factorize also 
the other correlations as follows: 



{V{t)V{t)) = 

(vmit)) = 
{vmm = 

It follows that 



(e^«*)(y^)(e-^«*) + (e'«* y)(ye-'«*) 
(\/e'^*)(y e~*^*) + e*^* y)(e-'^*) 
(e^^* y)(e-*«* V) + (e'^*)(y e-'^* V) 



C(0,0) = C(0) (65) 
C(t,t) - Poit)C{0) + coitf (66) 
C(0,t) = C(t,0) = co(t)2-co(t)co(t) (67) 

Substituting into Eq. ([56|) we get the result 

1/2 



AE^,,^{t) - (l+co"(t))C(0) - co(t)" + 2co(t)co(t) 



(68) 



We note that for short times (t ^ to) we can use the 
approximations 



co(t) 
co(t) 
co{t) 



-Cit) 



(69) 
(70) 
(71) 
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which demonstrates the agreement with Eq. ((58)) . The nu- 
merical results in Fig. ([TU)) confirm the validity of Eq. (p5)) 
for the FM, and Eq.dM]) for the WM. We note that in 
the FM case the effect of recurrences is more pronounced, 
because they are better synchronized: all the out-in-out 
traffic goes exclusively through the initial state. Fig.lfTTj) 
establishes the a/2 ratio throughout the whole range of 
s values. One should be aware that for small s the (w^/s) 
in Eq. ([5ni) should be replaced by (uj^/s) — (w^/s), which 
takes into account the finite level spacing in the numeri- 
cal simulations. For very small s this goes to \og{gujc), as 
if s ~ 0. The finite level spacing effect clearly shows up in 
the numerics, and would not arise in the strict continuum 
limit. 



IX. SUMMARY AND DISCUSSION 

In this paper, we have considered a quantum mechan- 
ical system, prepared in a discrete state, that subse- 
quently decays into a non-Ohmic continuum of other 
states. Two different models that have the same spectral 
properties, but still different underlying dynamics have 
been analyzed and the respective results have been criti- 
cally compared. One model (FM) refiects integrable dy- 
namics while the other is an RMT model (WM) that cor- 
responds to chaotic dynamics. In both cases, a universal 
generalized Wigner time that governs the relaxation pro- 
cess has emerged, while the non-universality is reflected 
in the decay law: We find that for "non-Ohmic chaos" 
(WM case) the survival probability becomes a stretched 
exponential beyond the Wigner time scale, which is both 
surprising and interesting. This is contrasted with the 
"integrable" power-law decay that takes over in the long 
time limit (FM case) , and obviously is very different from 
the Ohmic exponential result. Only the standard case of 
Ohmic bandprofile is fully universal. 

We have also investigated the temporal behavior of the 



second moment of the spreading wavepacket. We have 
found that in the FM case the generalized Wigner time 
is reflected in the spreading process and not only in the 
survival probability, contrary to the naive linear response 
theory expectation. 

Non-Ohmic coupling to the continuum emerges in 
various frameworks in physics. Quantized chaotic sys- 
tems might exhibit non-Ohmic fluctuations due to semi- 
classically implied long time power-law correlations, and 
in any case the typical power spectrum is in general not 
like "white noise" (e.g. [T^ - [lq V Other examples [l^ 
appear in the context of a many-particle system, where 
the hierarchy of states and associated couplings, ranging 
from the single-particle levels to the exponentially dense 
spectrum of complicated many-particle excitations, can 
lead to a very structural non-Ohmic bandproflle describ- 
ing the residual interactions. These non-universal struc- 
tures of the bandprofile of the perturbation, lead to a 
highly non-linear decay of the survival probability. De- 
pending on the context, the survival probability could be 
also related to the study of dephasing, or indirectly to the 
study of quantum fidelity and irreversibility: the general- 
ized Wigner time is reinterpreted as the coherence time, 
in the same sense as in Landau's Fermi liquid theory. 

It is worth mentioning, that in a bosonic second quan- 
tized language the decay of the probability can be re- 
interpreted as the decay of the site occupation n. If the 
interaction between the bosons is neglected this reduction 
is exact and merely requires an appropriate dictionary. In 
the latter context each level becomes a bosonic site which 
is formally like an harmonic oscillator, and hence the ini- 
tially empty continuum is regarded as a zero tempera- 
ture bath. Consequently the decay problem is formally 
re-interpreted as a quantum dissipation problem with an 
Ohmic (s=l) or non-Ohmic (s^l) bath. The generalized 
Wigner time scale is associated with the damped motion 
of the generalized coordinate h. 



model, s 


LDOS 


Po{t) 


FM, WM s = 1 


Lorentzian 


Exponential decay 


WM, 0<s<2 


Semicircle core + tails 


Exponential-like decay 


FM, 0<s<2 


Singular core 4- tails 


Exponential-like followed by Power-law 


FM, s = 2 


1/uj 


Log(f) decay 


FM, WM, s > 2 


Delta core + tails 


No decay after transient 


s < 


Core + low weight tails 


Infrared dominated 



TABLE L The various results for the LDOS and for the survival probability at a glance. 



Table 1 summarizes the various results that we have 
obtained for the survival probability. We conclude this 
section with a somewhat technical discussion of the 
crossovers between the various s regimes. Note that 



the statements below are implied by inspection of FigH) 
The strictly universal Ohmic result holds for s = 1. The 
super-Ohmic universal regime is bounded from above by 
s = 2, but for finite ultraviolate cutoff uJc the effective 
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border is a bit lower. This means that the marginally 
universal s = 2 behavior prevails in a finite strip around 
s ~ 2. For s > 2, if the coupling to the continuum is small 
enough, the survival probability does not decay, except 
a short transient that can be described by FOPT. Sim- 
ilarly, the sub-Ohmic universal regime is bounded from 
below by s = 0, but for finite infrared cutoff the effective 
border is a bit higher. Below this s-border the FOPT 
tails of the LDOS become sub-dominant [as implied by 
the divergence of the first term in Ea. ([55)) ]. and the de- 
cay of the survival probability becomes infrared deter- 
mined: This means that the effective cutoff is not tUc, 
but some different ill-defined (model dependent) cutoff 
at much lower frequencies, that might be determined by 
the level spacing statistics. 



Appendix A: Green function formulation 

For convenience we take the energy reference as 
i?o = 0. The Resolvent is defined as 



1 



oj-n + io 



(Al) 



Substitution oiH ~ Hq + V, expansion to infinite order, 
and exact geometric summation, can be carried out in 
the FM case, leading to the following standard result: 



iG+H|o) = 



1 



uj - A{u;)+i{T{uj)/2) 



(A2) 



Using the identity Im[G+] = —tt6{ijJ — H) one realizes 
that the LDOS is given by the expression 

pH = -ilm[(0|G+H|0)] (A3) 

leading to Ea. (p6|) . The evolution operator is given as the 
FT of the resolvent, namely U{t) FT [-2Im[G'+(a;)]], 
hence the survival amplitude is 

co(i) = (0|[/(0|0) = FT [27:p{io)] (A4) 

in agreement with the elementary derivation in the text. 



Appendix B: Optional derivation of the survival 
amplitude formula 

In this appendix we give an optional direct derivation 
for the survival amplitude formula in the FM case with- 
out relaying on the theory of Green functions. We are 
interested in finding c„(i) and co(t), which are the am- 
plitudes to find the particle in the respective levels. The 
Schrodinger equation is 



■ dco 

.dCn 



— iuJnt 



Vb,nC„(t) 



V„,oCo(i) 



(Bl) 
(B2) 



where ujn = En — Eq . By integration over (|B2p we get 



c„{t) 



VnfiCo{t')dt' 



placing Eq.dBll into Eq.dBl]) we get 



dt 



C{t-t') co(t')dt' 



where 



C(r)-FT[C(L.)]-^|K,ope-'""^ 



(B3) 



(B4) 



(B5) 



We want to solve the survival amplitude equation us- 
ing a Laplace transform technique. For that purpose we 
define K{t) = Q{t)C{t), where 6(r) is the Heaviside 
step function. Then we rewrite the equations as 



(B6) 



where co(t) is zero for negative r and unknown for posi- 
tive T. The corresponding equation for the Fourier com- 
ponents is 



(B7) 



From here it follows that the survival amplitude can be 
written as an FT 



co(t) 



27r 



^iuj + K{lu) 



(B8) 



From the definition of K{t) and using the convolution 
theorem it foUws that 



K{uj) = ^r{u}) + iA{uj) 
Hence consistency with Eq. ([T7| is established. 



(B9) 



Appendix C: The energy spreading formula 

For the derivation of the energy spreading formula 
Eq. (|56|) it is convenient to regard the Hamiltonian as time 
dependent. 



n{\{t)) =Ho + X{t)V 



(CI) 



where X{t) is a time dependent parameter. We define 
generalized forces in the standard way 



(C2) 



Using the Hcisenberg picture, and the usual notation 
J-{t) = U{t)^^TU{t), wc have the following relation 



dn 

'dt 



-mm 



(C3) 
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Thus the change in the energy can be written as 

Hit) - mO) = - f \{t')F{t')dt' (C4) 
Jo 

Squaring and taking the expectation vahie with respect 
to the initial state we get 

[l^E^,At)?= f dt' f dt"X{t')X{t")C{t',t") (C5) 

Jo JQ 

where C{t' ,t") = {F{t')F{t")) is the auto-correlation 
function. The free evolution during the time interval [0, t] 



corresponds formally to a "rectangular pulse" 

\{t') = e(t') - Q{t' - t) (C6) 

where 0() is a Heaviside step function. Its time deriva- 
tive is A = 5{t') - 5{t' - t), leading to Eg. ([56]) . 
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FIG. 3: (color online) Cartoon that illustrates the structure 
of the energy distribution ^(a;). The probability that would 
be held by the tails is po ^ 1 if FOPT strictly applied. It 
becomes of order unity (say 50%) for a fully developed core- 
tail structure. The core has a semicircle line shape in the WM 
case, and its border 7 is determined self-consistently. The 
distribution is characterized by Po{t), and by AiScore, and by 
AiSeprd, as explained in the main text. 



FIG. 1: (color online) Schematic illustration of the coupling 
matrix Vnm for the FM {upper panel) and the WM {lower 
panel). 



'0 




w 1 




FIG. 2: (color online) Illustration of the energy level scheme, 
and the definition of co and lj' , as explained in the main text. 



FIG. 4; (color online) The different ui regions of the LDOS 
are plotted for 0<s<2, relating to the (richer) FM case. The 
universal cu regions where the LDOS is cutoff- free, are shaded. 
The implied time-domain decay is indicated (note that upon 
the identification t ~ l/o; one can regard the horizontal axis as 
time stretching from right to left). The lower and upper cut- 
offs ujg and iJc, are indicated by vertical solid thick lines. The 
border of the different A{lj) expressions written in Eq. (|34p 
are plotted in curved solid (blue) lines. The curves of 70 and 
Svr/foo are plotted as well. The 70 curve was slightly modified 
for the purpose of presentation. The parameters used in this 
plot are a;e=l and uJc=W^ and e—2. 
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FIG. 5: (color online) Upper panel: The FM analytical LDOS 
for s=1.01. The core power-law (Ea. (|39l) ') is indicated by a 
dashed red line. Middle panel: The survival probability Po{t) 
computed using numerical FT of the LDOS. The exponential- 
like decay (based on Eq. (|49|) ) and the power-law decay (based 
on Ea. (|50|l 'l are plotted in dashed red lines. The vertical 
dashed black line emphasis the fact that the cross-over time 
too is different from to . Lower panel: The same for various val- 
ues of s £ [1.1, 1.5]. For short times there is a good agreement 
with an exponnential-like decay (based on Ea. (|49|l '). which is 
plotted in dotted black lines. 



FIG. 6: (color online) The LDOS for the FM and for the WM 
via direct diagonalization of 1600 x 1600 matrices with s = 1.5 
and e = 1.44. The units are such that Q=l, and hence the 
bandwidth is uJc = b. In the FM case formally b=N/2. Upper 
panel: The log-log scale emphasizes the emergent universal- 
ity of the tails as the cutoff tOc is taken to infinity. Lower 
panel: The log-linear scale emphasizes the difi'erence in the 
non-universal core component: FM has a singular core, while 
WM has a smooth semicircle-like core. 
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FIG. 7: (color online) The survival probability Po{t) for the 
FM (top) and for the WM (bottom), as a function of t/to, 
in log-log scale, for various values of e. The predicted FM 
power-law is presented as a dashed line. For all curves g — 1 
and s = 1.5. 
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FIG. 9: (color online) Upper panel: Examples for the time 
evolution of AiJcom for s=1.5 and 6=800 in the WM case. 
Lower panel: The extracted departure time versus the ex- 
tracted inverse saturation value. This scatter diagram demon- 
strates the validity of one parameter scaling. 
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FIG. 8: (color online) Upper panel: A re-plot of Fig. ([7]) using 
Y — — ln[Po{t)]/t and X = t in a, log-log scale, for represen- 
tative runs, showing that the decay in the WM case is de- 
scribed by a stretched exponential. The red bold dashed line 
has zero slope, corresponding to simple exponential decay for 
s=l. Middle panel: Additional curves for various values of 
s are plotted. The highlighted (thickened) segments demon- 
strate exponential-like decay. Lower panel: The power a is 
extracted by fitting to the form Po{t) — exp [— of the high- 
lighted segments of the middle panel. The expected power 
a=2—s is indicated by a dashed red line. 



FIG. 10: (color online) Scaled spread AS,p,d/(aj^e^/s)^/^ ver- 
sus scaled time ujct for the FM and the WM. The theoretical 
predictions Ea. (|24p and Ea. (|28|l are plotted for comparison. 
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FIG. 11: (color online) The scaled saturation value of AEsj,id 
for the FM (lower curves approach 1) and for the WM (upper 
curves approach V2). The saturation value falls down as s 
becomes smaller due to the finite level spacing. To establish 
the latter statement we re-plot the same data but re-scale the 
spread as AE,p,d/[{uj',e^ / s) - (cj^eVs)!'''^- 



